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The von Karman-Howarth equation implies an infinity of invariants corresponding to an infinity of different asymp- 
totic behaviours of the double and triple velocity correlation functions at infinite separations. Given an asymptotic 
behaviour at infinity for which the Birkhoff-Saffman invariant is not infinite, there are either none, or only one or only 
two finite invariants. If there are two, one of them is the Loitsyansky invariant and the decay of large eddies cannot 
be self-similar. We examine the consequences of this infinity of invariants on a particular family of exact solutions of 
the von Karman-Howarth equation. 



I. INTRODUCTION 



Results from recent laboratory experiments^ suggest that classes of homogeneous turbulence decay exist which 
are at odds with classical theory2,. As the general theory of homogeneous turbulence decay is based on invariants 
of the von Karman-Howarth equation^i^, these recent experiments call for a fresh study of what is true about these 
invariants. The present letter provides such a study in the context of decaying homogeneous isotropic turbulence. 
However, the assumption of isotropy could be dropped by following, for example, the method of Nie & Tanvees^. 



II. INVARIANTS OF THE VON KARMAN-HOWARTH EQUATION 



Starting from the von Karman-Howarth equation for decaying homogeneous isotropic turbulence^il, we show that 
it is possible to derive an infinite number of different invariants corresponding to an infinite number of different 
conditions at infinity. This equation is 



(1) 



where u' = u'(t) is the r.m.s. of the turbulent fluctuating velocity component u, u'^f{r,t) =< u{x,t)u(x + r,t) > and 
u'''k{r,t) =< {x,t)u{x + r,t) >, the brackets signifying an average over realisations or over the spatial coordinate x 
which is defined on the same axis as the velocity component u. Note that r > 0, that /(O, t) = 1 and that reflection 
invariance implies k{0,t) = 0. It is natural to assume that all derivatives of / and k with respect to r are not infinite 
at r = 0. 

Given suitable conditions at infinity, equation ([T|) can be used to calculate the rate of change of u'^ Ji^°° '"'"f^Tr rf'' 
for an infinite range of values of m and n. Repeated integrations by parts yield 
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where T^^p — 1 and T™_p = m(m — l)...(m — p) if p is a positive integer, and = m (note that m does not 
need to be an integer). These integrations by parts yield the right-hand side of ((2)| provided that m > n + 1, n > 0, 
lim,^oo(r'"-"fc) = and lim,_,oo(r™'"- V) = 0. The integral /+°° r^'^dr is finite if lim,^oo(r™-"+V) = 0. We 
make the assumption that f{r) and k{r) do not oscillate at infinity. 

Noting that r™_„_i - 4(m - ti - l)T^_„+i ^ {1 + n - m)[4r^'_„+i - T^_„] for all n>0, ^ simplifies to 
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By considering linear combinations of pairs of integrals J^^°° r"^ ^^dr and /q^°° ^ ^{ dr for which 
M, we can form an infinite number of invariants. There are two cases. One where M = m — n = 4, in which case 
there is no need to consider such linear combinations because [4TJ^_„^i — T^_„] = [4r™ — T™] — for any integer 
m > 4. This case immediately yields 

— u I r -TT^dr — (4) 
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for any integer n > under the conditions \imr^cxi{r^k) — and \imr^oc{r^ f) ~ 0. These conditions ensure that 
the quantity u'^ J^°" r^'^" ^-f^dr is both finite and independent of time for any integer n > 0. When n = 0, this 
quantity is the well-known Loitsyansky invarianliii^. Integrations by parts show that this quantity is proportional to 
the Loitsyansky invariant for any n > because lim,— »oo(»'^/) ~ and f{r) is assumed not to oscillate at infinity. 
The second case is for M 7^ 4. In this case the following linear combinations of integrals J^°° r^'^-frdr and 
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where m — n 



M 7^ 4, n and n' are non-negative integers such that n ^ n' and Cmuu' 
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under the conditions that M > 1, lim,._s.oo(r^^fc) — and limr_>oo(''^^^^/) ~ and that Imnn' is well-defined. Hence, 
the von Karman-Howarth equation admits an infinity of possible finite integral invariants depending on conditions 
at infinity. 

Whilst M does not have to be an integer, the smallest integer value of M for which such invariants exist is M = 2. 
The particular choice M = 2, n' = and n — I recovers the Birkhoff-Saffman invariant^i^ 



3/210 = u ' 
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The use of a single integral in this expression instead of the two integrals in equation ([5]) is significant because 
3r\f + r'^ = ^(rV) leads to 



3/210 = m'^ lim (r^/). 



(8) 



showing that /210 = if limr-*cx)('"^/) — 0, but also that /210 takes a finite value if defined as in (O rather than ([5]) 
and if lim, — >oo(f^/) is finite. 
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The Birkhoff-Saffman invariant ((7| can be generalised into an infinite series of invariants in two steps. Firstly, for 
any n > 1, define 



hno = y-' 

for which the following iterative relation holds 
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Hence, if f{r,t) ^ az{t){L{t) /rf (where L{t) is a length-scale and a-^L^ ^ 0) to leading order when r — 00, then the 
generalised Birkhoff-Saffman invariants /2nO are finite and their time-independence implies the time-independence of 
a-iL'^u'^ (and vice versa). As a second step, define 



l-Mnn' = ^ 
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for any M > 1 and any n 7^ n' . Noting that CMn(n'+i) ~ —{M + n' + l)CMnn' , one can derive the iterative relation 

, an' r 

lMn(n'+i) = -(M + n' + 1)Jm„„' + u'^ hm {r"+" +'^). (12) 

Under the same assumption that f{r,t) ~ a-i,{t){L(t) /r)^ to leading order when r — >■ cxa, (|lip and (|12p applied to 
M = 2 can now be used to show that all generalised Birkhoff-Saffman invariants /2nn' are finite and their time- 
independence is equivalent to the time-independence of a^L^u'^. 

Hence, our generalised Birkhoff-Saffman invariants lead to a conclusion previously reached by Birkhoff^ and 
Saffman^ on the basis of the constancy of ((Tjl alone. Namely, if /(r) ~ a3(L/r)^ as r — >■ 00, and if lim, — ioo{r^k) = 0, 
then 

^(a3L\'^) = 0. (13) 
dt 

As for any Af > 1 but different from 2 and 4, the case A/ = 3 corresponds to a new set of integral invariants. 
Similarly to the A/ = 2 case, we rewrite the invariants /3„o using only one integral, i.e. 
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The condition limr^cx>{r* f) — under which we established the constancy of /3„o implies /310 = 0. However, 
if /3n0 is defined as in (I14|l rather than ([5}, then it is permitted to relax this condition and assume instead that 
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f{r,t) a4{t){L{t)/rf (where L{t) is a length-scale and a4L'^ ^ 0) to leading order when r — > oo. In this case, and 
without forgetting the accompanying condition lim, — nx>{r^k) — 0, Isno is finite for all n > 1, and its invariance in 
time leads to 

^(a4L'u'2) = 0. (17) 

An eff'ectively identical argument to the one given above for l2nn' shows that all integral invariants I^nn' are in fact 
finite and time-independent under the conditions that f{r,t) ~ a4,{t){L{t) / r)'^ to leading order when r ^ oo and 
limr->cx)(?'^fc) ~ 0. Their time independence is also equivalent to (|17|l . 

The cases M > 4 are similar to the cases M = 2 and M — 3. In general, for any M > 1 such that M 7^ 4, we have 

(M + 1)7mio = lim (r*^+ V) (18) 



and 



hl{n + l)0 = -^'AfnO + Ca/ {„+ l)o1t'^ hni (t"""*"^^" ) • (19) 



For simplicity, we focus on iMno because the argument based on (lll|) and (|12p which we gave for M = 2 can be 
applied here to show that what holds for ImuO also holds for iMnn' ■ A re-definition of Jji/no in terms of a single 
integral (|ll|l instead of ([5]) allows the possibility for non-zero invariants of order M . Specifically, with such a re- 
definition, it is possible to assume ,f{r,t) ~ aM+i{t){L{t)/r)^'^^ (where L{t) is a length-scale and OM+ii*^^^ ^ 0) 
to leading order when r — >■ oo. In terms of the energy spectrum E{k) in Fourier space, this assumption takes the form 
E{k) ~ aA/+iM'^i(KL)" in the limit k -> when M > 1 (because^ E{k) = ^ Jf^ dr{3f(r) + rdf /dr)Kr sin Kr). 
Under this assumption and the accompanying condition limr^oo(r*^fc) = 0, /MnO is finite for all n > 1 and its 
invariance leads to 

'^;{aM+iL'"+\'^) = 0. (20) 



This proves a more precise version of the principle of permanence of large eddies given in p. 113 of the 1995 book 
by Frisch**. Note that (|20|) has already been obtained by Rotta" and Lundgren^'' by direct inspection but without 
noticing the integral invariants pip and therefore without the resulting systematic approach given here. 

We stress that M does not need to be an integer for equations ((6|, (|18p . (|19p and (|20ll to hold. However, as Rotta^ 
remarked, E{k) results from an integral over a spherical shell in wavenumber space^ so that any Af < 2 would imply 
that the spectral tensor^ (that is the Fourier transform of the velocity correlation tensor Rij =< 'Ui{x)uj{x + r) >) 
diverges as A; — >■ 0. We therefore limit the remainder of this letter to AI > 2. There is no a priori upper limit to M 
as the results of this section are valid for any M > 1. 
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III. CONSEQUENCES OF THESE INVARIANTS 



It is clear that we have an infinity of possible invariants depending on the asymptotic behaviours of /(r, t) and k{r, t) 
at infinity. Some of these invariants can also be expressed in terms of the velocity correlation tensor Rij , specifically 
in terms of its trace Ru which is a function of only r — |r| because of homogeneity and isotropy. In Batchelor's book 
on turbulence^, one can find the identity Rii{r) = u'^(3/ + r^) for homogeneous isotropic turbulence. Using this 
identity, one obtains 

//"OC poc 
r""^i?iidr = 47r / r'^'^ Rii{r)dr ^ 4tv{M - 2)u'^ r"/(r)dr + 4WMim (r*^+ V)- (21) 
Jo Jo 

for any M > 2. 

As noted by BirkhofS and Saffman^, this integral equals A-ku'"^ limr->oo when M = 2 and is finite if this limit 
is also finite. If this limit vanishes, then so does J Rudr, but in both cases J Rudr is an invariant. 

For any M > 2, J r^'~^ Rudr diverges in the case where \imr^cx,{r^^'^^ f) is finite but equals A-k{M — 
2)u'^ r^^ f{r)dr in the case where lim, — ioo{r^^'^^ f) — 0. Hence, with the exception of M — 2 and M = 4, 
J r^'~^ Riidr is not in general invariant, even though there are invariants iMnn' for every value of M > 2. (The case 
M = 4 corresponds to J r^Radv which is, in fact, the Loitsyansky invariant in a different guise.) 

We now show that, for conditions at infinity which are such that the Birkhoff-Saffman invariant is not infinite, 
either none or only one or only two invariants are finite. Assuming that there exists a number Mf > 2 for which 
limr-!.oo(''*^-'^''^^/) ~ ^Mf+iL^^f^^ ^ and a number Mg for which \imr^ao{r^^ k) — for any AI in the interval 
2 < M < Mg but limr_>cx)(r*^fc) / for any AI > Mg, then the following five possibilities present themselves for 
iMnn' redefined in terms of a single integral pif) instead of ([5]). 

(i) Mg < Mf and Mg < 4, in which case all invariants Imnn' = for AI < Mg and all Imuu' for AI > Mg are not 
invariant . 

(ii) Mg < Mf and Mg > 4, in which case all invariants iMnn' = for AI < Alg except the Loitsyansky invariant 
which is the single non-vanishing invariant, and all Imu for M > Alg are not invariant. In this case u'^ Io'°° r*' f{r)dr 
is the only non-vanishing invariant. 

(iii) 4 > Mg > Alf, in which case all invariants iMnn' ~ for M < Mf but invariant iM/nn' 7^ and all integrals 
^Mnn' with M > Mf diverge. 

(iv) Mg > 4 > Alf in which case all integrals Imuu' for which M < Mf are invariant but vanish and iMjun' 7^ and 
is invariant. 

In cases (iii) and (iv), aMj+iL'^'i'^^u'^ is the only non-vanishing invariant and Alf < 4. 

(v) Aig > Alf > 4 in which case there are only two non-vanishing invariants when Mf > 4, the Loitsyansky invariant 

6 



and iMfnn' , i-e. u'^ Jo'°° r'^f{r)dr and aMf+iL^'f^^u'^. When M/ = 4, u'^ r'^f{r)dr is the sole non-vanishing 
invariant. 

All in all, depending on conditions at infinity, either no finite invariants exist, or, if such exists, then either 
aMj.+iL'^'f'^^u'^ is the sole finite invariant with Mf < 4, or u'^ Jo'°° ''"^f{r)dr is the sole finite invariant, or 
u'^ J^°° r'^f{r)dr and aMf+iL^'f'^^u''^ (with Mf > 4) are the only two finite invariants. 

We close this letter by testing this conclusion on George-type self-preserving solutionsiii^^ of ([1} because of 
the recent claim that it might be possible to engineer self-preserving decaying homogeneous isotropic turbulence in 
the wind tunnel^. These solutions are of the form f{r,t) = f[r/l{t)] and k{r,t) = b{iy,u'(,,lo,t — to)K[r /l{t)] where 
u'q = u'{to) and lo = l{to). Introducing these forms into ((TJ, one obtains the solvability conditions ^it'^ ~ —2avu'^/l^, 
■^l^ — cu and 6 — I3v/{u'l) where a > 0, c > and /3 are numerical constants. It follows that 

-2a/c 

(22) 
and 

l^{t) = ll+cv{t-to). (23) 

If the conditions at infinity are such that no finite invariant exists, then no obvious constraint can be imposed 
on the exponent 2a /c and the rate of turbulence decay. However, in the case where the sole finite invariant is the 
Loitsyansky integral, then 2a/c — 5/2. In the case where the sole finite invariant is aMj+iL'^'i'^^u''^ with 2 < Mf < 4, 
then we can take L{t) — l(t) and the self-preserving form of / implies that aMj+i must be constant in time. We 
therefore get 2a/c = [Mf + l)/2 which lies between 3/2 and 5/2. 

Finally, when the conditions at infinity are such that u'^ r'^f{r)dr and aMj+\L^'i'^^u''^ (where Mf > 4) 
are both finite and invariant, then no George- type self-preserving solution of ([1]) is allowed because of the time- 
independence of aM^+i implied by such solutions. Noting that the contribution to J^°° r'^f{r)dr coming from small 
values of r is negligible, this conclusion is valid more broadly for any form of /(r) which is permissible by and the 
incompressible Navier-Stokes equations and which conforms with self-similar decay of large eddies^, i.e. for which 
f{r,t) ~ f[r/l{t)] if r is large enough and om^+i is time-independent as a result. Hence, if f{r) decays faster than 
as r — s- oo (i.e. E{k) drops faster than fc'* as fc ^ 0), and if the asympotic behaviour of the triple velocity correlation 
function is such that two finite invariants exist at once (case (v) above), then the decay of the large eddies cannot be 
self-similar. 



U (t) = Uq 



l+-r^(t-to) 
'0 
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IV. CONCLUSIONS 



A summary of main conclusions is in the abstract. The nature of turbulence decay depends critically on the 
asymptotic behaviour of the double and triple velocity correlation functions at infinite separations. There are four 
cases depending on whether Mf/Mg is larger or smaller than 1 and whether min{Mf,Mg) is larger or smaller than 
4. 

When Mf/Mg is larger than 1 and min{Mf, Mg) is smaller than 4 there are no finite invariants. When Mf/Mg is 
larger than 1 but min{Mf, Mg) is larger than 4 there is only one finite invariant and this is the Loitsyansky invariant. 
When Mf/Mg is smaller than f, there is either one or two finite invariants dependending on whether min{Mf,Mg) 
is smaller or larger than 4. In both cases aMj+iL^'^^^u'"^ is finite and invariant but when min{Mf , Mg) is larger 
than 4, Loitsyansky's u'^ r*'f{r)dr is a finite invariant too. 

Self-preserving turbulence decays in accordance with (1221) and (|23|) and the infinity of possible invariants permitted 
by IT} cannot determine the exponent in (|22p without prior knowledge of correlations between points in the turbulence 
which are extremely far apart. In fact, these correlations can even be such that no conclusion whatsoever can be made 
on the value of the exponent in (I22[) . and the relatively high values reported for this exponent in some wind tunnel 
experiments-^ cannot be ruled out theoretically without prior knowledge of these correlations. The self-preserving 
decay which seems to have been observed in some instances of fractal-generaged homogeneous turbulence^ suggests 
that Mf and Mg cannot be such that 4 < Mf < Mg in such instances of turbulence if it is isotropic. Research with 
many fundamentally different ways of generating turbulence^ needs to be carried out so as to gain some understanding 
of what determines conditions at infinity and whether they are all physically possible. 
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